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Abstract 

We introduce the Mathematica package SummerTime for arbitrary-precision com¬ 
putation of sums appearing in the results of DRA method. So far these results 
include the following families of the integrals: 3-loop onshell massless vertices, 
3-loop onshell mass operator type integrals, 4-loop QED-type tadpoles, 4-loop 
massless propagators [l|, 2,(3, 4J. The package can be used for high-precision nu¬ 
merical computation of the expansion coefficients of the integrals from the above 
families around arbitrary space-time dimension. In addition, this package can 
also be used for calculation of multiple zeta values, harmonic polylogarithms and 
other transcendental numbers expressed in terms of nested sums with factorized 
summand. 


1. Introduction 

Multiloop corrections within the Standard model and in general quantum 
field theory are of great interest now. In particular, this interest is connected 
with continuing search of New Physics both in collider experiments, such as 
LHC, and in high-precision spectroscopy measurement. 

Probably, the most effective approach to the multiloop calculations is the 
differential equations technique. It was introduced long ago (EL 0, 3 13 and has 
an impressive record of achievements. However, the differential equations do not 
work for one-scale integrals as their dependence on the scale is totally determined 
by dimensional considerations. On the other hand, the one-scale integrals are 
ubiquitous: they appear both in physical observables, like lepton anomalous 
magnetic moments or asymptotic R(s) ratio, and as boundary conditions in 
differential equations for multiloop integrals. 

Some time ago the DRA method of the calculation of the multiloop integrals 
has been introduced in Ref. [i]. This method is based on using the dimensional 
recurrence relations and analytical properties of the integrals as functions of 
the space-time dimensionality d. In a short time, the DRA method has been 
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successfully applied to a wide range of physically interesting one-scale families 

of integrals EL i EH, EH, 0, i, EH, EH, E3- 

The results of the DRA method are exact in d and have a form of nested 
sums with factorized dependence of the summand on the summation variables. 
Meanwhile, the applications require series expansion of the integrals in d near 
integer point, usually near d = 4. The most effective approach to obtaining 
analytical form of this expansion in terms of conventional transcendental num¬ 
bers proved to be the one based on high-precision calculation of the expansion 
coefficients and subsequent use of the PSLQ algorithm. Many analytical re¬ 
sults for the e-expansions of the master integrals obtained in this way have been 

already published in the papers EL i EH, EH S i, EH, EH, 


14] 


near 


d = 4 (in 

Ref. |10( also near d = 3). However, the results of DRA method contain much 
more information than the published expansion coefficients. In fact, they are 
quite analogous to the representation in terms of hypergeometric functions, and 
totally determine expansion coefficients around any value of d. The only prob¬ 
lem with these results is that, similar to the expansion of the hypergeometric 
functions, it may be quite difficult to express these coefficients in terms of con¬ 
ventional transcendental numbers, like multiple zeta values. Meanwhile, from 
time to time, there appears a necessity to calculate expansions of the master 
integrals either in different dimensionality or up to a higher order in e. 

The main goal of the present paper is to introduce a Mathematica pack¬ 
age SummerTime for arbitrary-precision calculation of the expansion coefficients 
of DRA results around arbitrary value of d. This package gives high-energy 
community the full access to the results of the DRA method. In addition, the 
package contains procedures for high-precision evaluation of the multiple-zeta 
values and harmonic polylogarithms. 


2. Solution of dimensional recurrence relation 

Dimensional recurrence relations seem to appear for the first time in Ref. EH 
and then rediscovered independently and systematically used in Ref. 16]. For 
the column-vector J of the master integrals of a given sector (i.e., the integrals 
with a given set of denominators) the dimensional recurrence relation has the 
form 

J(i/ + l) = C(i/)J(i/)+fl(i/), (1) 

where inhomogeneous term R(v) contains master integrals from the subsectors, 
v = d/2. Let S(z') be a revertible matrix, satisfying the equation 


S(i/) =S(i/+l)C(i/). 

Then the general solution of Eq. m can be written as follows: 

J( v) = § _1 (^) E§E)C ~ l (v)R(v). 

Here £ is the indefinite sum symbol with the property 
(E E)(zz + 1) — (£ F)(y) = f(v). 


( 2 ) 


(3) 


(4) 
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If the function F{v) decreases faster than 1/V when v —> +oo and/or v —> — oo, 
one can write 


+ OO 

E F(v) = E> F{y) + w(z) = f - 55 F(v + n) + u>(z) (5) 

n=0 

and/or 

-l 

T,F(u) = E< F(y) + uj(z) = f 55 F(y + n) + ui(z), (6) 

n =—oo 

where w(z) = tL>(exp(27rzz/)) is arbitrary periodic function. The key feature of 
the DRA method is the determination of the function cj(z) from the analytical 
properties on the integrals as functions of v. 

Assuming that simpler master integrals in R(v) are also represented in a 
form ® , we naturally come to the problem of calculation of the nested sums of 
the form 

I = s^ 1 Sr K Sk • • • S/ -1 Et 2 S 2 T 21 S^ S T1 SiTio§o 1 • ( 7 ) 

Here are rational rectangular matrices, t* £ {^,<}, and §i(v) is a 

hypergeometric matrix term, i.e. §~ (z^Si^+l) is a rational matrix. Explicitly, 
we have 

J(^) = 55 + n K) ■ ■ • (§2T 2 i§r 1 )(^ + n K + ■ ■ ■ + n 2 ) 

riiTiO 

x (SiTioSo 1 )(i / + tik + ... + ni) ■ (8) 

Here rZjTjO denotes rii < 0 or rii ^ 0 depending on whether r 7 ; is equal to “<’ or 
Passing to the variables ki = rii + TH+i + ■ ■ ■ + tik, we obtain 

!(«/) = §K 1 (v)§K(v + kK).-.(§2T21§T 1 )(v + k2)(.§lT 10 §o 1 Kv + kl). 

kiTiki -|_i 

( 9 ) 

In many physically interesting families of integrals each sector contains at 
most one master integral. Then all matrices are replaced by scalars and we 
come to the problem of evaluating the sums of the following form 

Jty) = 55 F K {y-\-k K )---Fi[y-\-ki), (10) 

kiTiki +1 


where each Fi is a hypergeometric term, i.e. the ratio Fi(v + 1 )/Fi{y) is a 
rational function, and kx+i = 0. In particular such are all the families of 
integrals calculated within DRA method 01 Hill Si El El El except for 
Ref. El where the result for multimaster integral was obtained. The present 
version of Summer Time package is restricted to the calculation of the scalar sums 

m- 
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3. Tree sums 


Remarkable property of the sums © and (in is that the dependence on 
the summation variables in the summand is factorized. If the limits of the sums 
were decoupled, the sum would be a product of one-fold sums. Nevertheless, 
even with the limits given by inequalities ktTiki+i, the sums in Eqs. © and 
m can be organized without nested loops as we explain below. But first we 
introduce the notion of tree sums and explain how sums of the form © or m 
can be rewritten in terms of those. 

Let us define the following correspondence between a certain vertex-labeled 
directed graph and the solution of a set of inequalities. Namely, let us consider a 
directed graph with nodes marked by some expressions. Then the corresponding 
set of inequalities contains inequality m ^ n for each edge m —> n. E.g., the set 
of integer triples (ki,k 2 ,k 3 ) satisfying inequalities (£3 < 0 )&(fc 2 ^ k 3 )k,(k\ < 
k 2 ) corresponds to the graph @^{^3 + l)«{fc 2 + l)-*(fc| + 2} All sums appearing 
in the DRA method can be labeled by path graphs. However, convergence 
rate analysis is easier for the sums which summation limits determined by the 
directed rooted tree graphs. Any path graph can be reduced to a set of directed 
tree graphs. E.g., 


[fc3 + 1J [fci + 2] 



which corresponds to the identity 


y, F (k 3 , k 2 , h) = y F(-n 3 - 1, n 2 , -n 1 - 1) 

&3 <0 0^713 

k 2 ^k3 0^n 2 

k\<k 2 O^rii 

+ y F(—n 3 - l,n 2 + l,m) -I- y F(—n 3 - 1, — n 2 - 1, — m - 2). (12) 

0^713 0^712 

O^Tli 712^^1 

711^712 712^713 

Note that if the summand on the left-hand side of Eq. CE2D was a hypergeometric 
term with factorized dependence on the summation variables, the same is also 
true for the summands on right-hand side of the equation. The factorized form 
of the summand is a very important property which allows one to organize the 
calculation of the sums without nested loops. 
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In what follows we will use term ‘tree sum’ for the sums with limits deter¬ 
mined by the directed rooted tree and with factorized dependence of the sum¬ 
mand on the summation variables. It is very instructive to use the following 
alternative recursive definition of tree sums: 

Definition 3.1. A tree sum is 

• Case 1: Any one-fold sum T(k ) = X]^Lfc a ( n ) 

• Case 2: Any sum of the form T(k) = Y^=k a ( n ) 11^=1 ^(n), where T,;(n) 
are tree sums (j is a natural number). 


3.1. Computation of tree sums 

The numerical estimate of the tree sum can be obtained by replacing in 
the above definition the upper limit oo by some sufficiently large number N. 
A naive prescription would give the computational complexity 0(N K ) of such 
an estimate, where K is a ‘nestedness’ of the sum (the number of summation 
variables). In order to organize the computation much more effectively, we 
use recursive nature of the definition l3Jl We first note that the evaluation 
of the list of one-fold sums {T(0),T(l),... ,T(N)} can be done in one path 
with the computational complexity O(N). Then, assuming that we already 
have such lists for each T), we can evaluate T(k) from Case 2 of definition 13.11 
also in 0(N). Therefore, the overall computational complexity of obtaining 
the estimate is O(NK). We stress the difference between the computational 
complexities of the sums with factorized and nonfactorized summand. 


Convergence acceleration. Once we eliminated the exponential dependence on 
K , we have to examine also the dependence of the complexity on the required 
precision P (P is a number of decimal digits). For exponentially converging 
sums (in what follows we call them G-sums ) we have N ~ P and the overall 
computational complexity it0 O(PK). However, for power-like behavior of the 
summand (we call such sums the H-sums) the number of terms in the naive com¬ 
putation depends exponentially on P. Therefore, we need to devise a suitable 
algorithm for convergence acceleration. 

Let us consider a one-fold sum 

:x ^ 

T = ^ a(n ), a(n) ~ —— at large n , a > 0 . (13) 

n =0 


Then in a naive prescription the required number of terms is N r^i 10 p/a . For¬ 
tunately, it appears that we can use the convergence acceleration described in 


Ref. [17j. The general idea of the convergence acceleration is to replace the 


1 We do not take into account the penalty of dealing with the P-digit numbers, which makes 

the complexity likely to be 0(P 1 2 K). 
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original sum with the modified one converging to the same value but with faster 
convergence rate. Let the transformed sum have the form 

OO 

T 1 = J2 a i(n) (14) 

n —0 

and 

ai(n) = (1 - f(n))a(n) + f(n + 1 )a(n + 1), (15) 

where f(n) is some function to be fixed and a(—1) = 0 by definition. Then the 
partial sum is expressed via as 

N 

T (N) = a\(n) = TW - f(0)a(0) + f(N + 1 )a(N + 1). (16) 

n =0 

The partial sums an TW both have the same limit when 

/(0) = 0, /(n)o(n)"^°0. (17) 

The convergence rate is boosted when 

ai(n)/a(n) " 1 - f{n) + f{n + 1)(1 + l/n ) _a_1 0 . (18) 

All these conditions are satisfied by f(n) = n/a. Note that if we know N terms 
of the original sequence a(n), we can reconstruct N — 1 terms of the transformed 
sequence. 

The leading asymptotics of the transformed term «i(n) depends on the next- 
to-leading term of a(n). Let us assume that the asymptotic expansion of a(n) 
has the form 

OO 

a ( n )~£j7Tfc' ( 19 ) 

k =1 

Then the transformed sequence has the expansion 

OO 

( 2 °) 

fe=i 

This allows one to apply the same acceleration technique recursively. If we have 
N terms of the sequence, we may apply the above acceleration technique N — 1 
times which gives an estimate N\gN ~ P. Taking into account that com¬ 
putational complexity of the acceleration is 0(N 2 ), we have a rough estimate 
0(P 2 / \g 2 P) < 0(P 2 ) for the overall complexity of the power-like decaying 
sums. Note that the convergence acceleration is compatible with our algo¬ 
rithm of evaluation of the nested tree sums. The only complication comes from 
the fact that convergence acceleration leads to strong numerical compensations 
when calculating the terms of the accelerated sequence, therefore the working 
precision W should be taken higher then the required one, and the empirical 
rule is W = 7 P, where 7 ~ 1.7. Then, for A'-fold sum we have 0(Kj K P 2 ) 
complexity estimate. 
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Expansion in e. Tree sums appearing in the DRA method depend on space- 
time dimensionality d. One is usually interested in their Laurent expansion 
near some point (typically, near d = 4). For a generic sum, depending on e, the 
Laurent expansion does not always commute with the summation. However, the 
e dependence of the tree sums appearing in the DRA method is very special: 
any summation variable n enters the summand only via combination n + e or 
n — e. Therefore, the sums converge uniformly in e which means that the Laurent 
expansion and the summation commute. Therefore, in order to calculate the 
former, one can first expand the summand and then make the summation. It is 
important that our convergence acceleration procedure is linear. Therefore we 
can apply it to formal power series in exactly the same way as to numbers. 


4. Conventional transcendental numbers and functions as tree sums 

The main purpose of the SummerTime package is to calculate sums which 
arise in the DRA method. However, it is remarkable that many conventional 
transcendental numbers and functions can be expressed and effectively evaluated 
as tree sums. Let us consider the Goncharov polylogarithms defined via 


Lir, 


,(®i ,---,Xk) = 


E 


nr* J- nr* * nr* 

•V J 2 • • • it/jl. 




Xr* . . . Xi 


( 21 ) 


To represent the above definition in the form of a tree sum, it suffices to rewrite 
it as 


ij+j 

x j 


Llni,...,n fe (a:i ,...,x k ) - E II (j . I j) 

Multiple zeta values 


( 22 ) 


C(n fc ,...,m)= n 

ifc>...>ii>0 j—1 


(sgnnj) lj 

„-KI 


are trivially related to the Goncharov poly logarithms: 

C (n k , ■ ■ •, ni) = Li| ni | i ...j„ fc |(sgn?7 i,..., sgn n k ), 


(23) 


(24) 


and, therefore can be readily represented in the form of tree sums. Harmonic 
polylogarithms (HPL) are defined by (see, e.g., Refs. 0,0) 


JT(0,...,0;aO = 


(lr 


n\ 


(25) 


and 


H(n,n k ,...,n i,x)= f„(t)H(n k ,... ,nnx), (26) 

J o 
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where n £ { 0 , ± 1 } and 


= /»(*) = -, = ( 27 ) 

It is standard to use nonzero integer n in the indices of I? as a shortcut for se¬ 
quence of \n\ — 1 zeros appended by sgnn, e.g., H( 3, —2; a;) = H( 0, 0,1, 0, —1; x). 
Since trailing zeros can be eliminated with the help of functional relations, see 
fl9l |. we may consider only H(n k , ■ ■ ■, n\] x), where all n\, ... ,n k are nonzero 
integers. We have 

H{n k ,..., nr, x) = (cti <j 2 ---Ok) Li|n 1 |,...,|n fe | (<7102, cr 2 cr 3 ,..., a k x), (28) 
where Oi = sgn rij. Therefore, HPLs also can be expressed via tree sums. 


5. Application example 


So far, the DRA method was applied to the following families of the integrals: 

• To 3-loop onshell massless vertices [lj; 

• To 3-loop onshell mass operator type integrals 

• To 4-loop QED-type tadpoles [lO]; 

• To 4-loop massless propagators ii- 


Partly to 4-loop onshell mass operator type integrals 14]. 


The results of the DRA method for the first four families can be downloaded 
as . zip files from the package web page. Together with the SummerTime package, 
they allow one to calculate with arbitrary precision the expansion of the integrals 
near any value of d. 

Let us present a minimal example of the application of SummerTime package. 
We assume that the content of the archive M. zip is put to a working directory of 
the Mathematica session. Then the following program calculates the expansion 
of M 6j 2 integral near d = 3 up to e 3 (e = (3 — d)/ 2) with the precision 1000 
digits: 


inti]:- <<SummerTime 1 

SetDirectory [NotebookDirectory[]] 

TriangleSumsSeries [<<"M62"/.d -> 3-2e,{e,3},1000] 
w[1] , - 126.330936333943790321081484798414734532015352412682120... x e “ 2 
+ 449.863309410101651375361946647344871810340547434784147...x e _1 
+ 7492.35397292225323538373824972066060688984257832532489815... 

- 17762.4348499407483506220786909230708553370113393194388... x e 

- 148136.431403928982879168608861673380547613003055965761... x e 2 





- 3048084.17207428077785717945708196781177068038009440576... x e 3 

+ 0(e 4 ) 

On the notebook with 2 cores Intel core i7 the run time of the above program 
was 12 min. 

In order to showcase the usability of the package, we have calculated the 
expansions of all four-loop massless propagator integrals near d = 3 (the cor¬ 
responding results for d = 4 can be found in Refs. HiS). The results are 
presented in Appendix. The calculation provides a strong evidence of appear¬ 
ance of new transcendental numbers in addition to conventional ones: multiple 
zeta values and alternating harmonic sums. 

6. Brief manual of Summer Time package 

The Summer Time package can be downloaded from the site 
http://www.inp.nsk.su/~lee/programs/SummerTime/ and installed accord¬ 
ing to the instruction included in a distribution archive. After the installation, 
the package can be loaded into Mathematica session by a command: 

<< SummerTime' 

The full list of functions provided by the package can be seen by ?SummerTime ' *. 
Here we will consider only the basic ones. 

6.1. Summation specifications 

The SummerTime package supports two types of summation limit specifica¬ 
tions. The first one, called triangle, is designed to handle the sums of form GOD. 
This summation specification represents a list of variables with ± signs, e.g.: 

{±n,±k,...}. (29) 

Signs determine the direction of summation: {..., ±n, k,... } denotes ±n < k < 
oo and {..., ±n, —k, ... } denotes n < k < oo. 

The second specification is used for tree sums, given by definition 13.11 It 
represents a tree, which is written in the following form: 

{node, childi, child 2 , ... }, (30) 

where children are also trees. If node doesn’t have any children, it’s simply 
written as: 

{node}. (31) 

The node itself is a summation variable. The root of a tree is always an integer. 
For example, the specification for a tree sum with limits 0 < n < k < oo can be 
written as { 0 , {n, {fc}}}. 
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6.2. Summation functions 

The SummerTime package defines following eight summation functions: 

TreeSum[expr, specs, p] 

TreeSumSeries [expr, specs, {e, o}, p] 

TriangleSum[expr, specs, p] 

TriangleSumSeries[expr, specs, fe, o}, p] 

TreeSums[exprs, p] 

TreeSumsSeries [exprs, {e, o}, p] 

TriangleSums [exprs, p] 

TriangleSumsSeries[exprs, {e, o}, p] 

Here expr denotes a summand, specs is a summation specification (either tree 
or triangle), exprs is a list of pairs {expr, specs}, e is a small parameter e, 
o is an expansion order with respect to e, and p is required precision. Func¬ 
tions names are constructed in the following way: Tree (Triangle) means that 
function uses tree (triangle) summation specification, Series indicates that the 
expansion of sum in e is performed, and Sums instead of Sum designates that sev¬ 
eral sums are calculated and their total is returned as a result. The advantage 
of Sums functions is that they perform computations in parallel. 

All summation functions have an optional boolean argument 
Progresslndicator, which denotes whether a progress bar must be shown dur¬ 
ing computations. Default value is True. 

Sum functions also have an optional argument called Information. By this 
parameter, user can set up a string identifier, which will be shown on top of the 
progress bar. This option can be helpful for a sequence of calls of the summation 
functions: using it, user can understand which particular sum is currently being 
calculated. 

Parallel functions (Sums) also have an optional boolean argument Parallel, 
which toggles, whether computation must be parallelized. Default value is True. 

6.3. Other functions 

Currently, the SummerTime package contains also implementation of multiple- 
zeta values and harmonic and Goncharov poly logarithms. By default these func¬ 
tions are disabled and in order to use them, one needs to set up a value of special 
variable $0ptions before loading the package: 

SummerTime‘Package‘$0ptions = {"MZV" -> True, "HPL" -> True, 

"MPL" -> True}; 

<< SummerTime' 

Multiple zeta values, harmonic and Goncharov polylogarithms are available via 

MZV[ml, ..., mk] £ (mi,..., m*,) 

HPL [{ml, ..., mk}, x] H (mi,..., m*,; x) 

MPL [{ml ,xl}, ..., {mk,xk}, x] Li mi) ... )TOfc (aq,..., x k ; x) 

Numerical values of these functions can be obtained via standard Mathematica 
function N. 


10 



7. Conclusion 


In the present paper we have introduced a Mathematica package SummerTime 
for calculation of the multiloop integrals obtained within DRA method. The 
main purpose of the package is an arbitrary-precision calculation of the ex¬ 
pansion coefficients of the integrals around any value of d. In addition, the 
package contains convenience tools for the calculation of various transcenden¬ 
tal constants and functions: multiple zeta value, harmonic polylogarithms and 
Goncharov poly logarithms. 

In the future, we plan to extend the package to the computation of the 
matrix sums. The difficulties connected with the convergence acceleration of 
the matrix sums seem to be possible to overcome using the approach of Ref. 
21]. Such an extension of the package is highly desirable given the families of 
integrals of actual interest for the present moment. 


Acknowledgments. This work was supported by RFBR grant No. 15-02-07893. 


Appendix A. Expansion of massless propagator-type master inte¬ 
grals 


Here we present the results for the expansions of the massless four-loop 
propagator-type master integrals around d = 3 (e = (3 — d)/2). It turns out [22| 
that these expansions may be of some interest for the applications. 

We fix normalization by dividing all expansions with 


Mqi (3 — 2e) 


r(i- e ) 5 r(4e-i) 

r(f-5e) 


(A.l) 


All results have been obtained by numerically calculating the expansions of the 
results of Ref. Q and then using Mathematica built-in function 
FindlntegerNullVector with a set including non-alternating and alternating 
harmonic sums. For some coefficients this approach failed which strongly indi¬ 
cates appearance of new constants in the expansions. We added the unresolved 
coefficients to the transcendental basis and used the extended basis for subse¬ 
quent integrals. In the end we have expressed all expansions in terms of standard 
transcendental numbers 


Cn = E 


1 




oo 


E 


i 

2 k k n ’ 


s 6 


EE 


(-i) 


m-\-k 


m 5 k 


(A.2) 


and the following ten unresolved coefficients: 

M 2 ( P = -262.1990984105955060153722678925879686999467800030..., (A.3) 
Mjf = 2873.2506195273103871316630429729277291187753911438..., (A.4) 
Af 2 ( f = -20850.77252988625214437845006030730940289184143944..., (A.5) 
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Af 2 (4) = 123806.28243923947828962358530545969884163429504554..., (A. 6 ) 

= 1915.2828601394850124013600200661647497224250709008 ..., (A.7) 
Af 2 (3) = -9230.990439724705220633706650782234978383599409697..., (A. 8 ) 
M 2 ( 2 } = 24179.184087947556033170242057981283994228407821870 ..., (A.9) 
= 20971.632987103828923594052378653898688768972653284..., (A.10) 
Aijf = 30385.287535356282840631913693445372174956734072988 ..., (A.ll) 
Afj 4) = -157029.840084536053310520481290249759485412480424.... (A.12) 


( n \ 

The notation of these constants is as follows: M A denotes the expansion coef¬ 
ficient of Ma( 3 — 2e)/.Moi(3 — 2e) in front of e n . 

We present results only for the integrals which can not be expressed in 
terms of T-functions. Where it is possible, we pull out rational factor making 
the expansion uniform in transcendental weight. 


M 21 ( 3-2e) 
Moi(3-2e) 


Af 22 (3-2e) 

Moi(3-2c) 


= + Af 2 ( 2 ) e 2 + M^e 3 + Af 2 ( 4 ) e 4 + O (e 5 ) 

= M$e + M 2 ( 2 ) e 2 + M^e 3 + M^e 4 + O (e 5 ) 


(A.13) 
(A. 14) 


M 23 (3 - 2e) 

M 0 i(3-2e) 


7 r 2 2 4 e ( 4 e - l)(10e - 3)(10e - 1) 
6e — 1 


1 + 20 £ 2 e 2 + 242^3 e 3 


+ 2487C 4 e 4 + (4840C 2 C3 + 18054Cs) e 5 + (29282C| + 231135Ce) e 6 + O (e 7 ) j 

(A.15) 


^ = Tr 4 2 8£ e (4 e -l)(10e-3)(10e-l)(l+48C 2 e 2 + 284C 3 e 3 +5706C 4 e 4 
Moi{6 — Ze) l 

+ (13632C2C3 + 19356Cs) e 5 + (40328C 2 + 496040Ce) e 6 + O (e 7 ) 1 (A.16) 


Af 25 (3-2e) 

Af 0 i(3-2e) 


7 r 2 2 4 e_ 1 ( 4 e — l)(10e — 3)(10e — 1) f _ 2 . 

-^--- l { 1 - + 74 &' 


- 213C 4 e 4 + (I 686 C 5 - 296C 2 C3)e 5 + (2738<f - 5367Ce) e 6 + O (e 7 ) \ (A.17) 


M 26 (3 - 2 e) 
Af 0 i(3-2e) 


(4e — l)(10e — 3)(10e — 1) 
6e+ 1 


|3C2 + (42C 3 -12 ai C 2 )e 
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( 32a 5 
32a?( 2 - 324aiC4 - —- 1 + 768a 5 + 42( 2 ( 3 

+ 1023( 5 ) e 3 + ( - 32a? ( 2 + 648a?( 4 - 168aiC 2 C3 + + 3072a 6 + 1086( 3 2 
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+ 


23967(6 


-2112s 6 ^e 4 + 0(e 5 ) | (A.18) 


M 27 ( 3-2e) (4e — l)(10e — 3)(10e — 1) 

M 0 i(3-2e) “ 6e- 1 


6( 2 + (48ai( 2 - 42( 3 )e + (96a 2 ( 2 
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